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ABSTRACT 

In the faint star KIC 9700322 observed by the Kepler satehite, 76 frequencies with 
amplitudes from 14 to 29000 ppm were detected. The two dominant frequencies at 9.79 
and 12.57 d~^ (113.3 and 145.5 /zHz), interpreted to be radial modes, are accompanied 
by a large number of combination frequencies. A small additional modulation with a 
0.16 d~^ frequency is also seen; this is interpreted to be the rotation frequency of the 
star. The corresponding prediction of slow rotation is confirmed by a spectrum from 
which 11 sin i = 19 ± lkms~^ is obtained. The analysis of the spectrum shows that 
the star is one of the coolest S Set variables. We also determine Toff = 6700 ± 100 K 
and log 3 = 3.7 ± 0.1, compatible with the observed frequencies of the radial modes. 
Normal solar abundances are found. An £ — 2 frequency quintuplet is also detected 
with a frequency separation consistent with predictions from the measured rotation 
rate. A remarkable result is the absence of additional independent frequencies down to 
an amplitude limit near 14 ppm, suggesting that the star is stable against most forms 
of nonradial pulsation. A low frequency peak at 2.7763 d~^ in KIC 9700322 is the 
frequency difference between the two dominant modes and is repeated over and over 
in various frequency combinations involving the two dominant modes. The relative 
phases of the combination frequencies show a strong correlation with frequency, but 
the physical significance of this result is not clear. 

Key words: stars: oscillations - S Set - stars: individual: KIC 9700322 - Kepler 



* Based on observations obtained with the Hobby-Eberly Tele- 
scope, which is a joint project of the University of Texas at Austin, the Pennsylvania State University, Stanford University, 
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1 INTRODUCTION 

The Kepler Mission is desig ned to detect Ea rth-like plan- 
ets around solar- type stars (jKoch et al.ll201ol ). To achieve 
that goal, Kepler is continuously monitoring the brightness 
of over 150 000 stars for at least 3.5 yr in a 105 square de- 
gree fixed field of view. Photometric results show that after 
one year of almost continuous observations, pulsation am- 
plitudes of 5 ppm are easily detected in the periodogram for 
stars brighter than V = 10 mag, while &t V = 14 mag the 
amplitude limit is about 30 ppm. Two modes of observation 
are available: long cadence (29.4-min exposures) and short- 
cadence (1-min expo sures). With short-cadence exposures 
l|Gilliland et al.ir201Gl ) it is possible to observe the whole fre- 
quency range seen in 5 Set stars. 

Many hundreds of 5 Set stars have now been detected 
in Kepler short-cadence observations. This is an extremely 
valuable homogeneous data set which allows for the ex- 
ploration of effects never seen from the ground. Ground- 
based observations of 5 Set stars have long indicated that the 
many observed frequencies, which typically span the range 
5 — 50 d~^, are mostly p modes driven by the K-mechanism 
operating in the Hell ionization zone. The closely-related 
7 Dor stars lie on the cool side of the 5 Set instability strip 
and have frequencies below about 5d~^. These are g modes 
driven by the convection-blocking mechanism. Several stars 
exhibit frequencies in both the 5 Set and 7 Dor ranges and 
are known as hybrids. iDupret et al.l |2003) have discussed 
how the K and convective blocking mechanisms can work 
together to drive the pulsations seen in the hybrids. 

The nice separation in frequencies between (5 Set and 
7 Dor stars disappears as the amplitude limit is lowered. 
Kepler observations have shown that frequencies in both 
the 5 Set and 7 Dor regions are pr esent in almost all of th e 
stars in the 5 Set instability strip (|Grig:ahcene et al]|2010l ). 
In other words, practically all stars in the 5 Set instability 
strip are hybrids when the photometric detection level is 
sufficiently low. 

Statistical analyses of several 5 Set stars observed from 
the ground have already shown that the photometrically ob- 
served frequencies are not distributed at random, but that 
the excited nonradial modes cluster around the frequencies 
of the radial modes over many radial orders. The observed 
regularities can b e partly explained by modes trapped in the 
stellar envelope (|Breger. Lenz fc Pamvatnvkhl 120091 ) . This 
leads to regularities in the observed frequency spectra, but 
not to exact equidistance. 

In examining the Kepler data for 5 Set stars we no- 
ticed several stars in which many exactly equally-spaced fre- 
quency components are present. There are natural explana- 
tions for nearly equally spaced frequency multiplets such as 
harmonics and non-linear combination frequencies. In some 
of these stars, however, these mechanisms do not explain the 
spacings. In these stars there is often more than one exact 
frequency spacing and these are interleaved in a way which 
so far defies any explanation. 

Some examples of equally-spaced frequency components 
which remain unexplained are known from ground-based 
observations. The 5 Set star 1 Mon has a frequency triplet 

Ludwig-Maximilians-Universitat Miinchen, and Gcorg-August- 
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where the departure from equidistance is extremely small: 
only 0.000079 ± 0.000001 d"^ (or 0.91 ± 0.01 nHz), yet the 
frequency splitting cannot be due to rotation because I = Q 
for the central component and 1=1 for the oth er two modes 
jBreger fc Kolenberell2006l : iBalona et aLlbOOll ) . In the P Cep 
star 12 Lac, there is a triplet with side peaks spaced by 
0.1558 and 0.1553 d~^. The probability that this is a chance 
occurrence is very small, yet photometric mode identifica- 
tion shows that two of these modes are 1 = 2 and the third 
is ^ = 1. This is therefore not a rotationally split triplet 
either (|Handler et al.ll2006l 'l. 

One solution to these puzzling equally-spaced frequen- 
cies coul d be non-linear mode intera c tion th rough frequency 
locking. iBuchler. Goupil fc HansenI (|l997l ) show that fre- 
quency locking within a rotationally split multiplet of a 
rapidly rotating star (150 to 200 kms~^) could yield equally- 
spaced frequency splitting, which is to be contrasted to the 
prediction of linear theory where strong departures from 
equal splitting are expected. 

In this paper we present a study of the 5 Set star 
KIC 9700322 (RA = 19:07:51, Dec = 46:29:12 J2000, Kp 
= 12.685). There are two modes with amplitudes exceed- 
ing 20000 ppm and several more larger than 1000 ppm. The 
equal frequency spacing is already evident in these large am- 
plitude modes. This star does not fall in the unexplained cat- 
egory discussed above. It is, however, a remarkable example 
of a star in which combination frequencies are dominant. 

The star has a large pulsational amplitude which 
can easily be observed from the ground. It was found 
to be variable in the "All Sky Automated Survey" 
(|Pigulski et all 120081 ). where it is given the designation 
ASAS 190751-1-4629.2. It is classified as a periodic variable 
(PER) with a frequency of 7.79 d"^ This is the 2 d~^ ahas 
of the main frequency (9.79 d^^), which is determined be- 
low from the Kepler data. The Kepler data is, of course, 
not affected by daily aliasing. It was also examined dur- 
ing t he "Northern Sky Variability Survey" (|Wozniak et al.l 
l2004h with up to two measurements per night. Due to the 
short periods of the star, the 109 points of NSVS 5575265 
were not suitable for a comparison with our results. 



2 NEW OBSERVATIONS OF KIC 9700322 

This star was observed with the Kepler satellite for 30.3 d 
during quarter 3 (BJD 245 5093.21 - 245 5123.56) with short 
cadence. An overview of the Kepler S c ience Processing 
Pipeline can be found in Ijenkins et al.l (|2010l ). The field 
crowding factor given in the KAC is 0.016, which is about 
the average for the Kepler field. The data were filtered by 
us for obvious outliers. After prewhitening the dominant 
modes, a number of additional points were rejected with a 
four-sigma filter as determined from the final multifrequency 
solution. 42990 out of 43103 data points could be used. We 
emphasize that most rejected points are extreme outliers 
and that the present conclusions do not change if no edit- 
ing is performed. As can be expected from near-continuous 
set of observations with one measurement per minute, the 
spectral window is very clean with the second highest peak 
at 0.046 d~^ and a height of 22% relative to the main peak. 

A small, typical sample of the Kepler measurements is 
shown in Fig.[T] Inspection of the whole light curve indicates 
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Figure 2. Two portions of observed spectrum matched to a model with T^g = 6700, log g = 3.7 (red line). In the left panel we show the 
region around H/3 that is sensitive to temperature, and in the right panel the region around Mgl triplet sensitive to gravity. 



that the pattern shown in Fig. [T] is repeated every 0.72 d. 
The repetition, however, is not perfect. This simple inspec- 
tion already suggests, but does not prove, that most of the 
variability is caused by a few dominant modes and that ad- 
ditional, more complex effects are also present. 



The (|Wozniak et all [20041 ') Input Catalogue also does 
not list any photometry for this star, but some informa- 
tion on the spectral energy distribution is available. The 
spectral energy distri bution was constructed using literature 
ph otometry: 2MASS (ISkrutskie et al]|2006l). GSC2.3 B an d 
R jLasker et al.ll2008l ). TASS V and I (jProege et"alll2006l ). 



and CMC14 r' (|Evans. Irwin fc Helmeill2002f ) magnitudes. 
Interstellar Na D lines present in the spectrum have equiva- 
lent widths of 60 ± 15 mA and 115 ± 20 mA for the Di and 
D2 hu es, respectively. The calibration of iMunari fc Zwitterl 
1^93) gives E{B - V)= 0.03 ± 0.01. 

The dereddened spectra l energy distrib ution was fit- 
ted using solar-composition (jKurucj Il993al ) model fluxes. 
The model fluxes were convolved with photometric fllter 
response functions. A weighted Levenberg-Marquardt non- 
linear least-squares fitting procedure was used to find the 
solution that mimimized the difference between the ob- 
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served and model fluxes. Since the surface gravity is poorly 
constrained by our spectral energy distribution, fits were 
performed for logg = 4.5 and logg = 2 to assess the 
uncertainty due to unconstrained \ogg. A final value of 
Tcff = 7140 ± 310 K was found. The uncertainties in TcS 
includes the formal least-squares error and that from the 
uncertainties in E{B — V) and log*;. We note here that in 
the next section values with considerably smaller uncertain- 
ties are determined from high-dispersion spectroscopy. 



3 CHARACTERIZATION OF THE STELLAR 
ATMOSPHERE 

In order to classify the star with higher precision and to test 
the very low rotational velocity predicted by our interpre- 
tation of the pulsation spectrum in later sections, a high- 
dispersion spectrum is needed. KIC 9700322 was observed 
on 2010 Au gust 12 with the High Resolution Spectrograph 
l|Tulllll998l ) on the Hobby- Eberly Telescope at McDonald 
Observatory. The spectrum was taken at i? ~ 30 000 using 
the 316g cross-disperser setting, spanning a wavelength re- 
gion from 4120 — 7850A. The exposure time was 1800 sees. A 
signal/noise ratio of 194 was found at 593.6 nm. We reduced 
the data using standard techniques with IRAF^ routines in 
the echelle package. These included overscan removal, bias 
subtraction, flat-fielding, order extraction, and wavelength 
calibration. The cosmi c ray effects were removed with the 
L.A. Cosmic package (jvan DokkumllioOll l. 

The effective temperature, Teff, and surface gravity, 
log g, can be obtained by minimizing the difference between 
the observed and synthetic spectra. We used a fit to the 
H/3 line to obtain an estimate of the effective temperature. 
For stars with Tes < 7000 K the Balmer lines are no 
longer sensitive to gravity, so we used the Mgl triplet at 
5167.321, 5172.684, and 5183.604 A for this purpose. The 
goodness-of-fit parameter, , is defined as 



N 



5 1 oh 



where A'^ is the total number of points and Jobs and 7th are 
the intensities of the observed and computed profiles, respec- 
tively. 5/obs is the photon noise. The error in a parameter 
was estimated by the variation required to change x'^ by 
unity. The projected rotational velocity and the microtur- 
bulence were determined by matching the metal lines in the 
range 5160 - 5200 A. 

From this procedure we obtained Toff = 6700 ± 100 
K, log,g = 3.7 ± 0.1, i;sini = 19 ± 1 kms'S ^ = 2.0 ± 0.5 
kms"\ In Fig. [H we show the match to observed spec- 
trum. The theoretical p rofiles were computed with SYNTHE 
Kurucz fc Avrettlll98ll ) using ATLAS9 atmospheric models 



Kurucj|l993bi r The solar opacity distribution function was 



used in these calculations. The effective temperature cal- 
culated from the spectrum is somewhat lower than that 



^ IRAF is distributed by the National Optieal Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 3. Abundance pattern derived for KIC 9700322. 



obtained by matching the spectral energy distribution dis- 
cussed in the previous section. The difference is within the 
statistical uncertainties. We note that the spectrum showed 
no evidence for the presence of a companion. 

Because of problems of line blending, we decided to use 
direct matching of rotationally-broadened synthetic spectra 
to the observations in order to determine the projected ro- 
tational velocity. For this purpose, we divided the spectrum 
into several 100 A segments. We derived the abundances in 
each segment using mi nimization. We used the line lists 
an d atomic paramete r s in Kurucz fc Belli (| 19951 ) as updated 
bv lCastefli fc Hubrid (|2004h . 

Table [1] shows the abundances expressed in the usual 
logarithmic form relative to the total number of atoms Atot. 
To more easily compare the chemical abundance pattern in 
KIC 9700322, Fig. [3] shows the steflar ab undances relative 
to the solar values (|Grevesse et ahllioiOl ) as a function of 
atomic number. The error in abundance for a particular el- 
ement which is shown in Table [T] is the standard error of 
the mean abundance computed from all the wavelength seg- 
ments. This analysis shows that the chemical abundance in 
KIC 9700322 is the normal solar abundance. 

The effective temperature determined for KIC 9700322 
makes the star one of the cooler 5 Set stars. Other pulsators 
with similar temperatu res are known, e.g., 67 00 K for p 
Pup jNetopil et al.ll2008l 'l and 6900 K for 44 Tan (|Lenz et al.l 
l2010h . 



4 FREQUENCY ANALYSIS 

The Kepler data of KIC 97 00322 were an alyzed with the sta- 
tistical package PERIQD04 (|Lenz fc Breger. 2005, '). This pack- 
age carries out multifrequency analyses with Fourier as well 
as least-squares algorithms and does not rely on the as- 
sumption of white noise. Previous comparisons of multi- 
frequency analyses of satellite data with other techniques 
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Figure 4. Amplitude spectra of KIC 9700322. The amplitudes (Tabic [2]l arc shown both linearly and logarithmically. The top panel 
illustrates the dominance of the two excited frequencies and that standard ground-based photometry would only detect these and a few 
combinations. 



Table 1. Abundances derived for KIC 9700322 expressed in term 
of logN{el)/Ntat 



c 


-3.6±0.3 


Sc 


-9.0 ±0.3 


Co 


-6.8 ±0.2 


o 


-3.1±0.2 


Ti 


-7.2 ±0.3 


Ni 


-5.8±0.3 


Na 


-5.8 ±0.3 


V 


-7.7±0.3 


Sr 


-8.6±0.1 


Mg 


-4.4±0.4 


Cr 


-6.5±0.2 


Y 


-9.4±0.2 


Si 


-4.6±0.1 


Mn 


-6.8±0.2 


Ba 


-9.2 ±0.4 


Ca 


-5.8±0.4 


Fe 


-4.6 ±0.2 







such as SIGSPEC (|ReeKeiJ20Q7|) have shown that PERI0D04 is 
more con servative in a ssigninR statistical significances, leads 
to fewer jPoretti et al.ll2009l ). and hopefully also fewer erro- 
neous, pulsation frequencies, but may consequently also miss 
some valid frequencies. 

We did not concern ourselves with small instrumental 
zero-point changes in the data since we have no method to 
separate these from intrinsic pulsation. Consequently, our 
solution contains several low frequencies in the region be- 
low ld~^ which may only be mathematical artefacts of in- 
strumental effects. The suspicion concerning the unreliable 
low frequencies is confirmed when comparing the present 
PERI0D04 results with those from other period search pro- 
grams and different data editing. 

Following the standard procedures for examining the 
peaks with PERI0D04, we have determined the amplitude 



signal/noise values for every promising peak in the ampli- 
tude spectrum and adopted a limit of S/N of 3.5. The value 
of 3.5 (rather than 4) could be adopted because most low 
peaks do not have random frequency values due to their ori- 
gin as combinations. This standard technique is modified for 
all our analyses of accurate satellite photometry: the noise is 
calculated from prewhitened data because of the huge range 
in amplitudes of three orders of magnitudes. 

After prewhitening 76 frequencies, the average resid- 
ual per point was 430 ppm. The large number of measure- 
ments (42990) lead to very low noise levels in the Fourier 
diagrams as computed by PERI0D04: 7ppm (0 — 10 d~^), 
4.7ppm (10 - 20d"^), 3.9ppm (20 - 40d"^), and 3.6 ppm 
(40 — 200 d~^). At low frequencies the assumption of white 
noise is not realistic. 

Our analysis was performed using intensity units (ppm). 
The analysis was repeated with the logarithmic units of mag- 
nitudes, which are commonly used in astronomy. The differ- 
ences in the results were, as expected, minor and have no 
astrophysical implications. The only small difference beyond 
the scaling factor of 1.0857 involved neighboring peaks with 
large intensity differences, in which the weaker peak was in 
an extended 'wing' of the dominant peak: the effects are 
numerical from the multiple-least-squares solutions. 

KIC 9700322 shows only six frequencies with ampli- 
tudes larger than 1000 ppm, of which only the two main 
frequencies are independent. Although a few ground-based 
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Table 2. Multifrequency solution of KIC 9700322 and identifications. Frequencies are given in cycles d ^ and also in /jHz. Amplitudes 
arc in parts per million (ppm). 



Frequency 


Amplitude 


Identification 


Comment 


Frequency 


Ampfitude 


Identification 




/.iHz 








d-i 










lb 0.00021 


± 0.002^ 


± 3 
















Main frequencies 


















9.7925 


113.339 


27266 


/i 


Dominant mode 












12.5688 


145.472 


29463 


/2 


Dominant mode 












0.1597 


1.848 


80 


f'i 


Rotation, causes 


combinations 










11.3163 


130.975 


27 


fi 


Quintuplet 












11.4561 


132.593 


145 


h 


Quintuplet 












11.5940 


134.190 


354 




Quintuplet 












11.7200 


135.648 


221 


h 


Quintuplet 












11.8593 


137.261 


112 


h 


Quintuplet 












Combination 


frequencies 


















n Qi o/i 
u.oiy4 


o.Dy / 


174 


2/3 




22.7723 


263.588 


15 


/4 


+ /5 


n A'7Ci^ 


O.040 


33 


3/3 




23.0501 


266.719 


39 


/5 


+ /6 


U.Dooo 


/ .oy4 


25 


4/3 




23.7187 


274.522 


27 




2/8 


u. / uyo 


S 91 1 


25 


/2 - /s 




24.0249 


278.066 


34 


/5 


+ /2 


U.y ( 4o 


IL.ZoZ 


27 


/2 - /e 




24.1628 


279.662 


35 


h 


+ /2 


111 97 


1 9 S7Q 
iZ.O ( y 


46 


/2-/5 




24.2888 


281.120 


52 


f7 


+ /2 


±.DDOD 


1 Q 9^4 


25 


/s - /l 




24.4282 


282.733 


57 


fs 


+ /2 


Z. / (Do 


Q9 1 


2633 


/2"/l 




24.9779 


289.096 


16 


2/2 


"/3 




33.981 


27 


/2 - /l + /3 




25.1376 


290.945 


2663 




2/2 


5.5526 


64.266 


89 


2/2 - 2/i 




27.9139 


323.078 


203 


3/2 


-/l 


7.0162 


81.206 


632 


2/1 - /2 




29.3776 


340.018 


191 




3/1 


9.4731 


109.643 


57 


/l - 2/3 




32.1538 


372.151 


479 


2/1 


+ /2 


9.6328 


111.491 


50 


/l -/3 




33.9554 


393.002 


15 


fe + fi 


+ /2 


9.9522 


115.188 


54 


/l +/3 




34.2207 


396.073 


16 


h + h 


4-/2 


10.1119 


117.036 


38 


/l + 2/3 




34.9301 


404.284 


536 






12.2494 


141.775 


34 


/2 - 2/3 




37.7064 


436.417 


329 




3/2 


12.4091 


143.624 


33 


/2-/3 




39.1701 


453.357 


16 




4/1 


12.7285 


147.321 


31 


/2+/3 




40.4827 


468.550 


34 


4/2 


-/i 


12.8882 


149.169 


22 


/2 + 2/3 




41.9464 


485.490 


23 


3/1 


4-/2 


15.3451 


177.605 


502 


2/2 - /l 




44.7227 


517.623 


114 


2/1 4 


~ 2/2 


16.8088 


194.546 


70 


3/1 - /2 




47.4989 


549.756 


81 


/i 4 


-3/2 


19.5850 


226.679 


2225 


2/1 




50.2752 


581.889 


83 




4/2 


21.2486 


245.933 


43 


/5+/1 




54.5152 


630.963 


35 


3/1 4 


-2/2 


21.3865 


247.529 


37 


/6+/1 




57.2915 


663.096 


58 


2/1 4 


-3/2 


21.5125 


248.987 


15 


/7 + /1 




60.0678 


695.229 


34 


/i4 


-4/2 


21.6519 


250.600 


20 


/8+/1 




62.8440 


727.362 


15 




5/2 


22.2016 


256.963 


26 


h + h~ h 




67.0840 


776.435 


19 


3/1 4 


-3/2 


22.3613 


258.812 


4902 


h + h 




69.8603 


808.568 


21 


2/1 4 


-4/2 


22.5210 


260.660 


19 


h + h + /3 














Other peaks 


in the amplitude spectrum 














0.0221 


0.256 


347 


fee 




13.2417 


153.260 


21 




/72 


0.0555 


0.642 


95 


fer 




14.6254 


169.275 


17 




/73 


0.1346 


1.558 


35 


fea 


/4-/5 and /5-/6? 


22.3250 


258.391 


20 




/r4 


0.3542 


4.100 


25 


/eg 




24.1464 


279.473 


30 




/75 


12.5347 


145.077 


82 


/70 




51.7521 


598.982 


14 




/76 


12.5837 


145.645 


30 


/71 















^ Accuracy of frequencies determined experimentally (see Section 4.1), independent of amplitude. The numbers apply only to 
unblended frequency peaks. Because of the high quality of the Kepler data, the frequency accuracy is much better than the resolution 
calculated from the length of a 30.3 d run. 
^ Determined by a multiple-frequency least-squares solution. 
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Figure 5. Echelle diagram of the detected frequencies using 
2.7763 d~^, wiiich is the difference between the two dominant pul- 
sation modes. All patterns can be easily explained through com- 
bination modes, rotational splitting and modulation (see text). 



be expressed in as very simple way through the equation 
/ = m/i ± n/2, where m and n are small integers. The 
fact that fi and /2 are the two modes with the highest 
amplitudes makes this approach also physically reasonable 
(see below). 

We also detect a frequency at 0.1597 d~^ (called /s). 
This frequency is important, since additional patterns are 
also seen: a number of peaks are separated by exactly the 
value of /a (see middle panel of Fig.[7| . 

Altogether, 57 frequencies can be identified as numer- 
ical combinations and multiples involving /i, /2 and fs by 
comparing the observed to the predicted frequencies. We 
can essentially rule out accidental agreements. Let us con- 
sider the combination frequencies at frequencies larger than 
3d~^, where the noise figures in the amplitude spectrum 
are reliable. For our identifications the average deviation 
between the observed and predicted frequency value is only 
0.00021 d~^. Such agreement is remarkable if one considers 
that the Kepler measurements used a time base of only 30 d 
and that 1/T — 0.03 d~^ where T is the time span between 
the last and first observation. The present result is typical 
for Kepler satellite data. 

If we use the least-squares frequency uncertainties cal- 
culated by PERIDD04, on average the observed agreement is 
44% better than predicted. However, such calculations as- 
sume white noise, w hich is not warranted. We can ad opt 
the formulae given in iKallinger. Reegen fc Weisj (|2008l ) for 
the upper limit of the frequency uncertainty to include 
frequency-dependent noise. We calculated signal/noise ra- 
tios in 5d~^ bins centred on each frequency with PERI0D04 
using the prewhitened spectrum. With this more realistic 
approach, the observed deviation of 0.00021 d~^ is exactly 
a factor of two lower than the statistical upper limit. This 
supports our identifications. 



campaigns lasting several years have succeeded in detect- 
ing statistically significant modes with smaller amplitudes, 
1000 ppm can be regarded as a good general limit. Observed 
with standard ground-based techniques, the star would show 
few frequencies. In all, we find 76 statistically significant fre- 
quencies. 



4.1 The observed frequency combinations 

Most of the detected frequencies can be identified as parts 
of regular patterns (see Fig.[7|. Visual inspection shows that 
the most obvious pattern is the exact spacing of 5/ — 
2.7763 d~^. This is confirmed by statistical analyses of all 
possible frequency differences present in the data. However, 
this pattern is not continued over the whole spectrum, but is 
present as different patterns, repeated and interleaved sev- 
eral times. Consequently, a simple explanation in terms of 
a Fourier series (e.g., of a nonsinusoidal light curve) is not 
applicable. 

Fig. [5] shows the Echelle diagram using 2.7763 d~^, 
which demonstrates the presence of remarkable patterns. In- 
vestigation of these patterns reveals that they originate in 
very simple frequency combinations and that the 2.7763 d^^ 
is only a marker of the true explanation: combinations of the 
two dominant modes /i and /2, as shown in Table[51 In fact, 
2.7763 d-^ = (/2 -/i). 

The frequencies shown in the top panel of Fig. [7] can 



4.2 The quintuplet 

Five almost equidistant frequencies in the 11 — 12 d~^ range 
are also present together with various combinations of these 
frequencies with /i and /2. This is shown in the bottom 
panel of Fig. [7] 



4.3 Explanation of the Echelle diagram 

We can now explain the patterns seen in the Echelle diagram 
(Fig.O in a simple manner. The vertical structures are the 
combination frequencies involving /i and/or /2. They are 
displaced from each other because different low integers of 
m and n (in the equation / — mfi ± n/2) are involved. The 
horizontal structures with a slight incline correspond to the 
frequencies separated by ~0.13 and 0.16d~'^, which are con- 
nected with the rotational frequency of 0.1597 d~^ through 
rotational splitting and modulation. The incline occurs be- 
cause the small frequency differences between adjacent fre- 
quency values must show up in both the x and y directions 
of the diagram. Details on the values and identifications of 
the individual frequencies are listed in Table|2] 
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Figure 6. Amplitude spectrum of KIC 9700322 after prcwhitcn- 
ing 76 frequencies. Note the low amplitudes of any additional 
pulsation modes present in the star. 



Table 3. The additional 1=2 quintuplet: Observed separations 



Frequency Separation from central frequency 
d-i d-i 



± 0.0002 


± 0.0003 


11.3163 


-0.2777 


11.4561 


-0.1379 


11.5940 





11.7200 


0.1260 


11.8593 


0.2653 


4.4 Additional frequencies 



A few additional peaks have been identified, which are not 
related to /i — /g in an obvious or unique manner. The lowest 
frequencies were already discussed earlier as probable zero- 
point drift and their values were dependent on how the data 
were reduced. The non-combination frequency at 51.75 d^^ 
has an amplitude of only 14 ppm. Calculation of the noise 
around the frequency gives a signal/noise ratio of 3.6, which 
makes it a very marginal detection. 

Fig. [6] shows the amplitude spectrum of the residuals 
after prewhitening of the 76 frequencies. No peak is statis- 
tically significant and the overall distribution of amplitudes 
is typical of noise. Nevertheless, we have examined the high- 
est (not significant) peaks, since a few of these may be real. 
Three peaks can be identified with expected values of addi- 
tional combination frequencies, e.g., a peak at 8.129 d~^ can 
be fit by 2/1-/5 at an amplitude signal/noise ratio of 3.0. 



5 DISCUSSION 

Although this star was selected because of its very clear 
exactly equal frequency spacing, it turns out that the fre- 
quency spacing is explained as simple combination frequen- 
cies arising from non-linearities of the oscillation. This is 
different from another class of 5 Set stars in the Kepler 
database which also show exact frequency spacings, but in 
a manner which is not at present understood. Examples of 
this strange class will be presented in a separate paper. 

What makes KIC 9700322 interesting is the remarkable 
way in which the large number of frequencies are related to 
the two main frequencies, /i and /2. This behaviour is very 



similar to the high amplitude 5 Set star KIC 9408694, also 
discovered in the Kepler database. The frequency patterns 
together with their amplitudes permit us to identify the dif- 
ferent frequencies and to provide physical interpretations. 

5.1 The dominant radial modes 

The period ratio of /i and /2 is 0.779. This is close to the ex- 
pected period ratio for fundamental and first overtone radial 
pulsation. The pulsation amplitudes of J Se t stars i ncreas e 
with decreasing rotation, e.g., see Fig. 5 of iBregeJ (|2000t ). 
Furthermore, high amplitudes occur mainly in slowly ro- 
tating, radial pulsators: in fact, the high-amplitude 5 Set 
(HADS) subgroup is defined on the basis of peak-to-peak 
amplitudes in excess of 0.3 mag. Nevertheless, a rigid sepa- 
ration between radial HADS and lower-amplitude nonradial 
5 Set stars does not exist. Dominant radial modes with am- 
plitudes smaller than 0.3 mag have previously been found . 
Examples o f EE Cam ((Breger. Rucinski fc ReegenI 120071 ') 
and 44 Tau (|Lenz et al.ll2010l ). The situation might be sum- 
marized as follows: Dominant radial modes occur only in 
slowly rotating stars. 

Since KIC 9700322 is sharp-lined {vsmi = 19kms"^) 
and presumably also a slow rotator, it follows this relation- 
ship. The presence of two dominant radial modes with ampli- 
tudes less than the HADS limits of peak-to-peak amplitudes 
of 0.3 mag is not unusual. 

The measured vsmi value supports the interpretation 
of the observed 0.1597 d~^ peak as the rotational frequency. 
In fact, both dominant modes have very weak side lobes with 
spacings of exactly the rotational frequency. The side lobes 
are very weak: for /i and /2 the amplitudes are only 0.0018 
and 0.0011 of the central peak amplitudes. We interpret this 
as a very small modulation of the amplitudes with rotation. 
An alternate explanation in terms of rotational splitting of 
nonradial modes is improbable because rotational splitting 
does not lead to exact frequency separation unless there is 
frequency locking due to resonance. Also, the extreme am- 
plitude ratios tend to favour the interpretation in terms of 
amplitude modulation. 

Based on this mode identification assumption we in- 
vestigated representative asteroseismic models of the star. 
We have used two independent numerical packages: the first 
package consisted of the current versions of the Warsaw-New 
Jerse y stellar evolution code and the Dziembowski pul sation 
code (|Dziembowskilll977l : [Diiembowski fc GooASm^ . The 
second package is composed by the evolutionary c ode CE- 
SAM (Morel 1997I'). and the oscillat ion code filou (|Suarej 
l2002l : ISuarez. Goupil fc Morelll200l ). Both pulsation codes 
consider second-order effects of rotation including near de- 
generacy effects. 

The period ratio between the first radial overtone and 
fundamental mode mainly depends on metallicity, rotation 
and stellar mass. Moreover, the radial period ratio also al- 
lows for inf e rence s on Rosseland mean opacities as shown in 
iLenz et al.l (|2010l ). 

Indeed, an attempt to reproduce the radial fundamen- 
tal and first overtone mode at the observed frequencies with 
the first modelling package revealed a strong dependence 
on the choice of the chemical composition and the OPA L 
vs. OP opacity data (|lglesias fc Rogerd[l99^ : ISeatonll20"05l '). 
The best model found in this investigation was obtained with 
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OP opacities and increased helium and metal abundances. 
Unfortunately, this model (Tefi = 7400 K, logL/L© = 1.27, 
logg = 3.87, 2M0) is much hotter than observations in- 
dicate. The disagreement in effective temperature indicates 
that this model is not correct despite the good fit of the 
radial modes. 

As an additional test, by ad opting the radial linear 
non adiabatic models dev eloped bv iMarconi fc Paiial (liggsD 
and |m 

arconi et all (|2004 ). we are able to reproduce the val- 
ues of the two dominant frequencies with pulsation in the 
fundamental and first overtone modes, but with a lower pe- 
riod ratio (0.770) than observed. The best fit solution ob- 
tained with these models, for an effective temperature con- 
sistent with the spectroscopic determination and assuming 
solar chemical composition, corresponds to: A'l — 1.65M0, 
logL/Lo = 1.1, Tcff = 6700 K, logg = 3.83. We no- 
tice that for this combination of stellar parameters, both 
the fundamental and the first overtone mode are unsta- 
ble in these models. Moreover, looking at the Main Se- 
quence and post-Main Sequence evolutionary tracks in the 
gravity v ersus effective tempera ture plane, as reported in 
Fig. 4 of ICatanzaro et all (|2010l '). the solution Tcff = 6700 
K, \ogg = 3.83 is consistent with a 1.65M0 stellar mass. 

However, as already noted, the period ratio in our mod- 
els is lower than the observed value. To resolve this discrep- 
ancy, the possibility of low metallicity and rotation effects 
was examined in more detail with the second modeling pack- 
age. Models between T^b = 6200 K and 8600 K with masses 
between 1.2 and 1.76M0, were found to represent a good fit 
of /i and /2 as radial fundamental and first overtone, re- 
spectively. The best fit with the observations was found for 
M = 1.2 M0 models computed with omlt = 0.5, dov = 0.1, 
and a metallicity of [Fe/H] =-0.5 dex. Such a low value for 
the convec tion effic iency is in good agreement with the pre- 
dictions bv lCasas^t al. ( 200(^ ) for 5 Set stars, based on their 
non-adiabatic asteroseismic analysis. All these parameters 
roughly match the general characteristics of the S Set stars 
with dominant radial modes and large amplitudes, despite 
being in the limit in metallicity. 

The Pi/Po period ratios predicted by these models 
(which simultaneously fit Po) are near 0.775, which is lower 
than the observed ratio, 0.779. A period ratio of 0.775 is 
also obtain ed by adoptin g the radial linear nonadiabatic 
models by iMarconi et"aLl (|2004l ) at Z = 0.006, accord- 
ing to which the best fit solution with effective tempera- 
ture consistent with the spectroscopic determination, cor- 
responds to ^ = 0.006, Y = 0.25, M = 1.5, logL/L0 = 
1.04, Te = 6700 K, logp = 3.83. Again the fundamen- 
tal and first overtone modes are predicted to be simulta- 
neoulsy unstable for this parameter combination. We ex- 
plored the possibility that such a discrepancy might be due 
to rotation effects, p articularly second-order di s tortio n ef- 
fects, as discussed by Suarez. Garrido fc Goupill (|2006l ') and 
ISuarez. Garrido fc Moval (l2007l ). These investigations ana- 
lyze theoretical Petersen Diagrams including rotation ef- 
fects (Rotational Petersen Diagrams, hereafter RPDs), and 
show that Pi/Po ratios increase as stellar surface rotation 
increases. The rotation rate derived from observations is 
slightly below 25kms~^ (see section [5^ . At such rotation 
rates near degeneracy effects on the period ratio are small 
(less than 0.001 in Pi/Pq). However, when non-spherically 
symmetric components of the centrifugal force are consid- 



ered, near-degeneracy effects may be larger, around 0.0025, 
causing the presence of wriggle s in the RPDs (see Fig. 5 in 
ISuarez. Garrido fc Moval |2003) and Fig. 6 in jPamvatnykhl 
( 2003h V Such effects are even more significant for rotational 
velocities higher than 40 - 50kms~^. Consequently, near- 
degeneracy effects may help to decrease the discrepancy be- 
tween the observed period ratio and the slightly lower values 
predicted by the models. 

If the star had a low metal abundance (close to Pop. H), 
a detailed analysis of RPDs might have provided an indepen- 
dent estimate of the true rotational velocity (and thereby of 
the angle of inclination). However, the spectroscopic analysis 
indicates that the star has a solar abundance. KIC 9700322 
therefore represents a challenge to asteroseismic modeling, 
since it appears impossible to reproduce all observables si- 
multaneously with standard models. 



5.2 The combination frequencies 

We have already shown that the 50+ detected frequency 
peaks can be explained by simple combinations of the two 
dominant modes and the rotational frequency. Several differ- 
ent nonlinear mechanisms may be responsible for generating 
combination frequencies between two independent frequen- 
cies, vi and 1/2 ■ For example, any non- linear transformation, 
such as the dependence of emergent flux variation on the 
temperature variation (L = crT'*) will lead to cross terms 
involving frequencies + i>2 and vi — V2 and other com- 
binations. The inability of the stellar medium to respond 
linearly to the pulsational wave is another example of this 
effect. Combination frequencies may also arise through res- 
onant mode coupling when vi and V2 are related in a simple 
numerical way such as 2v\ « ?,V2- 

The interest in the combination frequencies derives from 
the fact that their amplitudes and phases may allow indirect 
mode identification. For nonradial modes, some combination 
frequen cies are not allowed depending o n the parity of the 
modes jBuchler. Goupil fc HansenI Il9971 ) which could lead 
to useful constraints on mode identification. Since /i and 
/2 in KIC 9700322 are both presumably radial, there are no 
such constraints. 

The identification of /i and /2 with radial modes allows 
us to investigate the properties of the Fourier parameters 
of the combination modes with the aim to disentangle less 
obvious cases and/ or solutions with a smalle r num ber of 
combination terms. iBuchler. Goupil fc Hanseiil (|l997h show 
that a resonance of the type /c = ni/i -I- ng/g leads to a 
phase 0r = 0c — (?ii0i + ^2 02). In the same way we may 
define the amplitude ratios Ay = Ac/{AiA2). To investigate 
how (j)r and Ar behave with frequency, we first need the best 
estimate of the parent frequencies. We obtained these by 
non-linear minimization of a truncated Fourier fit involving 
/i , /2 and all combination frequencies up to the 4"^ order. 
The best values are /i = 9.792514 and /2 = 12.568811 d'^ 
The resulting amplitude and phases are shown in Table|4] 
together with the values of (jir and Ar- The phases were 
calculated relative to BJD 245 5108.3849 which corresponds 
to the midpoint of the observations. 

Fig. |8] shows how Ar and (f)r vary with frequency. From 
the figure we note that Ar is largest for /i + f2, 2/i, 2/2 
and /2 — /i and very small for the rest. It is also interesting 
that (j)r is a relatively smooth function of frequency, being 
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Figure 7. Top panel: Amplitude spectrum of those frequencies which are numerical combinations of the two dominant modes only. 
This demonstrates the richness of the combinations. Notice that the patterns are more complex than a simple Fourier sum. The middle 
panel demonstrates that some of those combination frequencies show further splittings with (rotation) . The bottom panel shows an 
additional set of five frequencies together with the combinations of the set with /i and /2 . 



practically zero in the vicinity of the parent frequencies, de- 
creasing towards smaller frequencies and increasing towards 
higher frequencies. This result is almost independent of the 
choice of /i and /2. T he standard deviation of fi and /2 
is 0.0001 d"^ using the iMontgomerv fc Q'Donoghud 
formula. One may arbitrarily adjust /i and /2 in opposite 
directions by this value, and using the corresponding cal- 
culated values of the combination frequencies, fit the data 
to obtain new phases. The resulting (f)r versus frequency re- 
mains monotonic, but the slope does change. The smooth 
monotonic nature of the (j)r versus frequency diagram re- 
mains even for a change of ten times the standard deviation 
in opposite directions for /i and /2 and for arbitrary changes 
in epoch of phase zero. The result is clearly robust to obser- 
vational errors, but it is not clear what physical conclusions 
may be derived from this result. The behaviour is certainly 
not random and must have a physical basis. Note that for 
simple trigonometric products, will always be zero. 

Finally, we note that the amplitudes of the combi- 
nation modes relative to the amplitudes of their parents 
can be compared wi th values detected in the star 44 Tan 
iBreger fc Lenj[2008l ). They agree to a factor of two or bet- 
ter, suggesting that KIC 9700322 is not unusual in this re- 
gard, just more accurately studied because of the Kepler 
data. 



5.3 The quintuplet 

In addition to the quintuplet structure around the two dom- 
inant modes another quintuplet with different properties is 
present in KIC 9700322 (see the listing of fi to /g in Ta- 
ble|3]). The average spacing between the frequencies in this 
quintuplet is slightly smaller than the rotational frequency 
(0.1338d"^ vs. 0.1597d"^). This makes this quintuplet dif- 
ferent from the quintuplet structures found around the two 
dominant modes, which exhibit a spacing that corresponds 
exactly to the rotation frequency. Moreover, the distribution 
of amplitudes within the third quintuplet is fundamentally 
different to the patterns around /i and /2. The given char- 
acteristics support an interpretation of the quintuplet as an 
I — 2 mode. 

The location of the quintuplet near the centre in be- 
tween the radial fundamental and first overtone mode rules 
out pure acoustic character. Consequently, the observed 
quintuplet consists of mixed modes with considerable ki- 
netic energy contribution from the gravity-mode cavity. For 
such modes theory predicts a smaller (and more symmet- 
rical) rotational splittings compared to acoustic modes due 
to different values of the Ledoux c onstant C„i- Using the 
fram ework of second order theory (|Dziembowski fc Goodd 
Il992h we determined the equatorial rotation rate which pro- 



Regularities in frequency spacings 11 



Table 4. Best fitting amplitudes, Ac (ppm), and piiases 4>c (radi- 
ans), for tiic parent frequencies and their combination frequencies 
up to fourth order. Relative phases, (pr and amplitudes, Ar are 
also shown. 
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Figure 8. The relative amplitudes, Ar, and phases, (pr (radi- 
ans) as a function of frequency (d~ ^ ) . The dotted lines show the 
location of the parent frequencies. 



vides the best fit of tlie observed quintuplet with an £ — 2 
multiplet. Tlie best results were obtained for an equatorial 
rotation rate of 23kms~^. This is only slightly higher than 
the observed v sin i value of 19 km s~ , and therefore indi- 
cates a near-equator-on-view. The Ledoux constant, C„;, of 
the i = 2 quintuplet is 0.164. For quadrupole modes C„; 
ranges between ~0.2 for pure gravity modes to smaller val- 
ues for acoustic modes. With (1 - C„;) = 0.836 this leads to 
a rotational frequency, Urot = of around 0.16 d~^. Con- 
sequently, this theoretical result confirms the interpretation 
of /s as a rotational feature and of the quintuplet as / = 2 



modes. Further support is provided by the fact that we see 
various combinations of the quintuplet with /i and /2. 

Moreover, the location of the quintuplet allows us to 
determine the extent of overshooting from the convective 
core. In the given model we obtained Oov ~ 0.13 but the 
uncertainties elaborated in Section [5.11 currently prevent an 
accurate determination. 



5.4 Further discussion 

A remarkable aspect of the star is the fact that so few pulsa- 
tion modes are excited with amplitudes of 10 ppm or larger. 

In the interior of an evolved S Set star, even high- 
frequency p modes behave like high-order g modes. The large 
number of spatial oscillations of these modes in the deep in- 
terior leads to severe radiative damping. As a result, nonra- 
dial modes are increasingly damped for more massive 5 Set 
stars, which explains why high-amplitude(5 Set stars pulsate 
in mostly radial modes and why in even more massive clas- 
sical Cepheids nonradial modes are no longer visible. 

In general, we do not expect the frequencies in the 
5 Set stars observed by Kepler to be regularly spaced be- 
cause, unlike ground-based photometry, the observed pul- 
sation modes are not limited to small spherical harmonic 
degree, I. For the very low amplitudes detected by Kepler 
we may expect to see a large number of small-amplitude 
modes with high I. The observed amplitudes decreases very 
slowly with I and, all things being equal, a large number 
of modes with high I might be expected to be s een in 
5 Set and other stars dBalona fc Dziembowskil Il999l ) . The 
^Sct stars HP 50844 (iPoretti et al.l |2009| ) and HD 174936 
(|Garci'a Hernandez et al.| 20091 ) observed by CoRoT show 
many hundreds of closely-spaced frequencies and may be ex- 
amples of high-degree modes. The relatively small number 
of independent frequencies detected in KIC 9700322 stands 
in strong contrast to the two stars observed by CoRoT. 

It should be noted that, unlike many 5 Set stars ob- 
served by Kepler, KIC 9700322 does not have any frequen- 
cies in the range normally seen in 7 Dor stars. The only 
strong frequencies in this range are a few combination fre- 
quencies. Although we have identified significant frequen- 
cies below 0.5 d~^, it is not possible at this stage to verify 
whether these are due to the star or instrumental artefacts. 
At present, we do not understand why low frequencies are 
present in so many 5 Set stars. 

Regularities in the frequency spacing due to com- 
bination modes have already been observed from the 
ground even in low a mplitude 5 Set s tars. An exam- 
ple is the star 44Tau jBreeer fc Lend I2OOS). Fig. 2 of 
fereger. Lenz fc Pamvatnvkhl (|2009l ') demonstrates that all 
the observed regularities outside the 5 — 13d~^ range are 
caused by combination modes. For combination modes the 
frequency spacing must be absolutely regular within the lim- 
its of measurability. This is found for KIC 9700322. 
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